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Abstract 



We present a few estimates of energy densities reached in heavy- 
ion collisions at the CERN SPS. The estimates are based on data and 
models of proton-nucleus and nucleus-nucleus interactions. In all of 
these estimates the maximum energy density in central Pb+Pb in- 
teractions is larger than the critical energy density e c ~ 0.7GeV/fm 3 
following from lattice gauge theory computations. In estimates which 
we consider as realistic the maximum energy density is about 2e c . In 
this way our analysis gives some support to claims that deconfined 
matter has been produced at the CERN SPS. Any definite statement 
requires a deeper understanding of formation times of partons and 
hadrons in nuclear collisions. We also compare our results with im- 
plicit energy estimates contained in earlier models of anomalous J/tp 
suppression in nuclear collisions. 



1 Introduction 



There exist plenty of models and Monte Carlo generators of proton-proton 
(pp), proton-nucleus (pA) and nucleus-nucleus (AB) interactions. These 
include models based on strings and their fragmentation, e.g. [D El El E| , on 
partonic cascades E] , on hadronic cascades IHj , on combined parton and 
hadron degrees of freedom j^, and on other pictures of the initial state of 
the collision. After fixing a few parameters, these models are able to find a 
reasonable agreement with data. 

Another class of models for pA interactions is based on successive col- 
lisions of the incident proton with those nucleons in A which are present 
within a tube in A given by the trajectory of the proton in the nucleus 
and by the non-diffractive cross-section for proton-nucleon interaction, e.g. 
[TUl HIl H3 Uni H31 HSl HZl CHI 113 I201- The model can be naturally 
extended to AB interactions. There it is based on the Glauber picture of 
colliding tubes of nucleons, see e.g. [HI EE21 EH d EEH1 EH HZ|- Models of 
this type proceed at every step in accordance with the data available from 
pp and pA interactions. A crucial input comes from data on the formation 
time of hadrons and on the Drell-Yan process [21j . Such a method keeps the 
parameters under control and gains some information about the space-time 
evolution of the process. 

One of the most important quantities of interest in heavy-ion collisions is 
the highest energy density reached. This quantity is relevant to the possible 
approach to the quark-gluon plasma. In this paper we estimate the energy 
density reached at the CERN SPS by using a simple version of the Glauber- 
type model with tube-on-tube interactions. 

We shall compare our estimates with recent lattice gauge theory results 
|22] 125] . The lattice results have brought a new and most interesting in- 
formation on the type and parameters of the phase transition between the 
quark-gluon plasma (QGP) and the hadron gas (HG). The phase transition 
is most likely of a cross-over type, with critical temperature T c pa 173 MeV 
and the rather low critical energy density of e c = e(T c ) ~ 0.7GeV/fm 3 . 

The paper is organized as follows. In Section [21 we describe a simple 
version of the Glauber-type model of AB interactions, which we use in our 
calculations of rapidity distribution of the energy AE/Ay. We also present 
our estimate of the volume occupied by the central rapidity region of Ay pa 1 
right after the tube-on-tube interaction is finished. The energy density is 
then estimated as the ratio of AE/Ay and AV/Ay. More details of the 
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Figure 1: Left: geometry of non-central nuclear collisions. Right: layout of 
tube-on-tube interaction (plotted without Lorentz contraction). 

model are explained in Section El In particular we discuss the formation 
time of hadrons and its relationship to the space-time picture of subsequent 
energy losses of incident nucleons in a tube-on-tube interaction. In Section 0] 
we present our results. Comments and concluding remarks are deferred to 
the last section. 

2 The model 

For the sake of simplicity we shall take the nuclei as hard spheres with radii 
Ra = 1.2A 1 / 3 fm and homogeneous number density p = 0.138 fm -3 . In the 
transverse plane, the impact parameter is denoted as b and a point in the 
transverse plane of the nucleus A is specified by the transverse coordinate s. 
The angle between b and s*is denoted as 9. The situation is sketched in Fig.^ 
In a Glauber model, the first part of the nuclear collision is described as a 
sum of tube-on-tube interactions (Fig.Q). In our model we will be interested 
in the energy density contained in such tubes just after the interaction. 

The cross-section of both tubes is equal to the non-diffractive nucleon- 
nucleon cross-section a = 30 mb = 3fm 2 . The lengths of the tubes and 
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2L B are given by 



2L A {s) = 2^/ R\ - 



(1) 



2L B (b,s,6) = 2^R|T 



b 2 - s 2 + 2bs cos6. 



(2) 



The average numbers of nucleons in both tubes are 



(n A (s)) = 2L A (s)pa, (n B (b, s, 9)) = 2L B (b, s, 9)pa . 



(3) 



For this type of models three ingredients have to be specified: 

(i) The probability distribution of the number of nucleons within the col- 
liding tubes Pa{ u a) and P B (n B ). 

(ii) For given ha, n B the nucleons in the tube in A can be numbered, start- 
ing with the head of the tube as i — 1, 2, . . . , ua, and similarly in B 
j = 1, 2, . . . , n B . It is assumed that every nucleon in the tube in A 
collides with every nucleon in the tube in B. We have to specify the ra- 
pidity distribution of all nucleons before and after every nucleon-nucleon 
collision. 

(hi) Finally, we have to specify the production of secondary particles and 
compute their energy distribution in rapidity after every nucleon-nucleon 
collisions. Consider the collision of the i-th nucleon in the tube in A, 
with the j-th nucleon in the tube in B, which we refer to as an 
collision. Both nucleons have lost a part of their rapidity in interactions 
prior to the collision. If we denote the incoming rapidities in such 
a collision by y\ and y B , we need to specify dN (y A , y 3 B ) / dy dp 2 - or at 
least the pr-integrated distribution dN(y A ,y B )/dy of secondaries pro- 
duced in the (i, j) collision. From the rapidity spectrum one computes 
the rapidity distribution of energy AE^ / Ay of the produced particles 
within the interval —0.5 < y < 0.5. 

After having specified the items (i)-(iii), the energy contained in all sec- 
ondary particles produced from a collision of two tubes in the rapidity interval 
—0.5 < y < 0.5 is obtained as 



AE sec 




Ay 



PA(n A )P B (n B )J2Yl 



Ay 



(4) 
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In order to obtain the total energy within the given rapidity interval we have 
to add the energy of incident nucleons in the two colliding tubes, which end 
up in the rapidity interval —0.5 < y < 0.5 when the tube-on-tube collision 
is finished. We denote this contribution by the index "st" from nucleon 
stopping and obtain 

AE tot _ AE sec AE st 
Ay ~ Ay + ~AY ' ^ 

This is the total energy resulting from a tube-on-tube collision. We will have 
to divide it by the volume occupied by quanta which were produced from 
these two tubes. 

We shall now describe the three inputs (i)-(iii) in our model. 

(i) Number distribution. The distribution of the number of nucleons in 
both tubes is assumed to be Poissonian, with mean values fi A = ( n A(s)) = 
2La{s)p<j and similarly for fi B 

p / \ {^aT a exp(-/i A ) {fi B ) nB exp{-p B ) 

Pa\P>a) = ; , PbsP'Bj = j • (6) 

n A \ n B \ 

(ii) Rapidity loss. In each nucleon-nucleon collision the rapidity loss of 
both nucleons is Ay. We will show results calculated for Ay = 0.5 and 
Ay = 0.7. In the CMS of nucleon-nucleon collisions at the CERN SPS, the 
absolute value of the rapidity of incident nucleons is y = \y A \ = \y%\ ~ 3. In 
every collision it decreases by Ay. When the rapidity of a nucleon becomes 
less than Ay, the nucleon does not participate in further collisions and its 
energy contributes to AE st in Eq. (jSJ). The term AE st / Ay is thus estimated 
as 

AE st °° 

— r — = P Ai n A) Pb(jib) [nA,siow(n A , n B ) + n BMow (n A , n B )\ x 1 GeV , 

n A ,n B =0 

(7) 

where riA,siow and n BtS i ow are numbers of incident nucleons which end up in 
the final state with —0.5 < y < 0.5. 

(iii) Particle production. We assume that the production of secondaries 
is the same as it is in vacuum and use the parametrization due to Wong and 
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Lu jT3] for the computation of the energy of secondary particles produced in 
the collision. In this model, the rapidity of charged particles produced 
in a collision of nucleons with rapidities y\ and y B is given as 



dy 

where 



A ((!-*+)(! -*_))" , (8) 



x+ = exp ^ ~ y%A ^ x ~ = ex v(y 3 B - y)> ( 9a ) 



a 



3.5 + 0.7 In ^57;, (9b) 



m^T = l m; + J B|) 1 /2 ) ( 9c ) 



,2 



A = 0.75 + 0.38 In v ^~, (9d) 
B T = 0.27 + 0.037 In (9e) 

In Eqs. (|9j) . m n is the pion mass and the nucleon mass. The average 
transverse momentum of produced pions Bt is taken in units of GeV/c and 
the CMS energy of the (i, j) collision in units of GeV. As discussed in 
|13j . the parameters in Q were tuned by comparison with experimental data 



The energy of neutral secondary particles (mostly pions) in the final state 
is taken into account by assuming that n^o = (n w + +n 7r -)/2 and multiplying 
the factor A in Eq. (j9d|) by 3/2. In order to go from particle distribution to 
energy distribution we multiply the right-hand side of Eq. (jHj) by the mean 
pion energy (E^) and obtain 

^ = 1.5A((l-x + )(l-X-)y(E 7T ). (10) 

Here (E x ) = (ml + B\ + p|) 1/2 , with p 2 L = (p^)/2 = B\j2 for pions in 
the rapidity interval —0.5 < y < 0.5. The value of AE^/Ay calculated by 
Eq. (fTU|) is then inserted into Eqs. (JIJ) and ©. 

Our interest here is in the energy density of quanta which form the sys- 
tem when the tube-on-tube interactions are finished. These quanta include 
secondary particles produced according to Eq. (jHJ) and slowed-down nucleons. 
The energy density is given by 



AE tot /Ay 
V(L A ,L B ,Ay) 



(11) 
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where V(La, L b , Ay) is the volume occupied by the central rapidity unit 
when the collision of the two tubes is finished. It is only the volume of parti- 
cles involved in a single tube-on-tube process, not the total fireball volume. 
We estimate it as 

V(L A , L B , Ay) = (2L A / 7 + 2L B /j + 2v t )a . (12) 

Here, 7 is the Lorentz contraction factor, for Pb+Pb collisions at the CERN 
SPS 7^9. The third term 2voto stands for the delay due to formation time. 
We assume that the formation of particles effectively sets in after the two 
colliding tubes have crossed each other. In our calculation we put vq = 0.5 
and varied the parameter t . 

3 Comments on the model 

The Glauber model of pA and AB interactions has been studied by many 
authors [m [UJ H21 [131 [HI CS1 CHI [HI CHI CH EE! [23 EH [SHI [2ni EDI- Experi- 
mental data have been reviewed by Busza and Ledoux [3T|. As pointed out in 
|3~T] , one of the problems of this field is caused by insufficient accuracy of the 
data. The situation can change soon, since the NA49 Collaboration at the 
CERN SPS has recently obtained new data on pA interactions, which were 
only briefly published so far [221 EH] • Note that these data are in the same 
energy region where anomalous J/ip suppression [S3] occurs and increased 
production of multistrange baryons [35] has been seen. For a more precise 
formulation of the model one would need accurate data on nucleon stopping 
in pA interactions, on the production of secondary particles in multiple col- 
lisions of a proton in the nucleus, and on the formation time of secondary 
particles. 

We shall now discuss the assumptions made in our simple model as well 
as the choice of parameters. 

First we turn to our assumption about the rapidity loss in every nucleon- 
nucleon collision. In pp collisions, the proton rapidity loss is as large as 
Ay « 1 (21 EH and this is assumed in many models [Po [ IT7 | I2K | IT31 I2%[ I3U] in 
which the data on pp interactions are directly extended to multiple collisions 
in pA and AB interactions. In [TSJ [TH] the rapidity loss is put equal to 
Ay w 0.5 and in some other models ^U] the rapidity loss per collision is as 
small as Ay « 0.3. Recent data of the NA49 Collaboration [33] support our 
choice of Ay = 0.5. In order to see the influence of a larger Ay, we have also 



6 



made calculations with Ay = 0.7 which is closer to the assumptions made in 
most models. 

A crucial assumption made in our model is the way in which the volume 
V(La, Lb, Ay) is determined, see Eq. (fT2|). We discuss this point in more 
detail here. 

An important contribution to the volume in Eq. comes from the 

formation time. We want to note that in a pA interaction the collision of the 
incident proton with the tube in the nucleus is a very complicated process 
which we are unable to describe in detail. A description of this process by 
multiple collisions only gives the final state but does not make statements 
about the real intermediate stages of the process itself. In early critical 
comments on this type of description [36j it was already pointed out that a 
simple classical explanation contradicts the data on Drell-Yan production in 
pA and AB interactions. The argument goes along these lines: suppose that 
in the first collision in the nucleus the proton loses some part of its momentum 
and the parton structure functions immediately adapt to this change. In such 
a situation the cross-section for the production of Drell-Yan pairs cannot be 
proportional to A a with a very close to 1. The same problem appears in the 
case of Drell-Yan production in nuclear collisions, where the cross-section is 
again proportional to (AB) a with a close to 1. This issue has been recently 
discussed by Gale, Jeon and Kapusta [2T] • They have introduced a coherence 
time of the energy loss of the proton in pA interactions. This coherence time 
gives the delay after which the proton energy is degraded to the value seen 
in the final state. By analysing the data [37, on Drell-Yan pair production 
in pA interactions, they concluded that the average proper coherence time 
to = 0.4 ± O.lfm/c. They have also pointed out that this coherence time 
is related to the formation time of secondary hadrons in pA interactions. 
Indeed, when a secondary final-state hadron is able to interact with other 
particles, its energy must be already felt as lost by the incident proton. This 
estimate of the formation time may rather be considered as a lower limit, 
since a part of the data in can be explained by shadowing corrections 
and by the energy loss of incident partons when traversing the nucleus |3*H] . 
Possible effects of shadowing and energy loss corrections in the Drell-Yan 
data from Fermilab were discussed for the first time in j^Hj. Recent analyses 
can be found in |3~Hj . 

Although the model of pA interactions uses hadronic degrees of freedom, 
this does not imply their dominance in the intermediate stages of the process. 
Since the model refers only to the final state, it is quite possible that the 
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dynamics of the intermediate stages is dominated by partonic degrees of 
freedom. The relationship between intermediate and final stage is thus given 
by some form of parton-hadron duality. 

Let us come back to our estimate of the volume. In our scenario, the 
process of tube-on-tube collision can be finished only after both tubes traverse 
each other completely. The effective volume is further increased by the delay 
needed for the formation of secondaries, which move with rapidities up to 
y = ±0.5. The value of the formation time to is assumed to be of the order of 
1 fm/c and is a free parameter. This estimate is based on works quoted in [21] , 
in particular on ^Uj- For experimental data, see [JT]. The data underlying 
these estimates are rather old and a new and more accurate information on 
the value of the formation time is most desirable. 

Before finishing this section, let us discuss some alternative estimates of 
the volume. One could argue that most of the secondaries are produced in the 
early nucleon-nucleon collisions and therefore the total lengths of the tubes 
2La and 2L B should be replaced by some effective lengths, say l e A s ~ La and 
If w L B . Then 

\/ half (L A , L B , Ay) = (La/ 7 + W7 + 2v *o)<t . (13) 

On the other hand, the process is not finished when only the front halves 
of the tubes cross each other, and the production of secondary particles is 
disturbed by nucleons that come later. 

A very popular Bjorken scenario 42J has been devised for asymptotic 
energies. In this regime the nuclei are really contracted to pancakes of van- 
ishing thickness. In the CMS of Pb+Pb interactions at the SPS the nuclei 
are contracted to 2R^/ r y m 1.58 fm and neglecting this length is not realistic. 

Both these estimates would lead to smaller volumes than formula ()12j) and 
hence to larger energy densities. In what follows, however, we only report on 
results obtained with the estimate (|12jh which we consider as most realistic. 

4 Results 

We calculated the energy density as a function of the length of the tubes. 
Before plotting our results it is useful to quote realistic values for the max- 
imum length of tubes that fit into a few selected nuclei. These lengths are 
given by 2L A iax = 2Ra = 2.4v4 1//3 and are summarized in Tabled 
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Table 1: Largest length of tubes for some selected nuclei in units of fm. 



A 


12 C 


16Q 


32g 


64 Cu 


115 In 


108 Ag 




208p] 3 238JJ 


<2J^max 


5.49 


6.05 


7.62 


9.6 


11.67 


11.43 


13.65 


14.22 14.87 



In all our calculations the following choice of parameters is made: a = 
30 mb, 7 = 9, p = 0.137 fm -3 , y° A = 3,y° B = -3. 

In Fig. El we plot six curves corresponding to the critical energy density 
e = e c = 0.7GeV/fm~ 3 as a function of 2La, 2L B - They were obtained for 
Ay = 0.5 and 0.7, and different values of the formation time parameter t - 

The influence of stronger nucleon stopping on the results is very weak: 
there is only little modification in the curves as seen from Fig. El This 
demonstrates that the major contribution to the energy density comes from 
produced particles, i.e. the first term in Eq. (0). Quite naturally, longer 
formation times lead to larger volumes and lower energy densities, and the 
critical energy density is thus reached in collisions of longer tubes. 

The results are very optimistic. According to Tabled critical energy den- 
sity is just reached in the centre of head-on S+S collisions and certainly in 
S+Pb or S+U interactions, because the combination of longest tube lengths 
for these pairs of nuclei falls into the region above the critical curve in Fig. |21 
One is tempted to claim the existence of the QGP. However, if we accept the 
conjecture that anomalous J/ip suppression is connected with plasma produc- 
tion, our statement comes out too optimistic. Anomalous J/ip suppression 
was only observed for larger collision systems |34j . 

J/i/j suppression as a signature of QGP formation in nuclear collisions 
has been proposed by Matsui and Satz [3H] more than 15 years ago. The 
anomalous J j ip suppression has been discovered by the NA50 collaboration in 
1996 [33]. About a year later phenomenological models of the QGP formation 
and J/ip suppression in nuclear collisions at the SPS have been proposed 
by two groups: Blaizot, Ollitrault 03] (we will refer to this model as BO) 
and Kharzeev, Lourengo, Nardi and Satz (KLNS, |45 ). As discussed in 
Refs. |4*4^ and by Nogova et al. 0H]> the condition for QGP formation in 
the BO model can be stated simply as 

[BO] 2L A (s) + 2L B (b, Sl 6) > 23.5 fm. (14) 
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Figure 2: Lines of constant energy density e = e c = 0.7GeV/fm 3 calculated 
for t = 0.5, 1.0, 1.5fm/c. Nucleon stopping was set to Ay = 0.5 (thick lines) 
and Ay = 0.7 (thin lines). Other parameters are given in the text. Energy 
density is calculated by Eqs. (fTTj) and (|T2"|) . 

For the KLNS model the condition for the formation of QGP reads 

[KLNS] 2L A x 2L B > 5.87 fm x (2L A + 2L B ) . (15) 

In these equations, 2L A , 2L B are given in units of fm. The critical curves 
for plasma formation according to these two models are shown in Fig. El 
Note that for interactions of equal-length tubes they both lead to almost 
the same results. Differences can only show up in non-central collisions and 
interactions of different size nuclei. 

We compare with these curves the results of our energy density estimate 
with reasonable parameters values. The critical curves of BO and KLNS 
models correspond to an energy density about 1.8 times above e c . We inter- 
pret this result as a consequence of delayed thermalization and rapid expan- 
sion in both the longitudinal and transverse [UJ direction which leads to fast 
cooling. In order to have a collectively behaving plasma that will be able to 
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Figure 3: Curves in the 2La, 2Lb plane, which indicate the limits of anoma- 
lous J/i/j suppression in the BO jUj and KLNS j3H] models, given by Eqs. (f"H|) 
and (|15|) . respectively. These curves are compared with constant energy den- 
sity curves with e = 1.7e c and e = 1.9e c , which were calculated in our model 
for Ay = 0.5 and to = 1 fm/c. 

screen the cc interaction, particle production must lead to an energy density 
higher than e c . Only in such a case can it remain in a deconfmed state until 
the collective behaviour is established. 

5 Comments and conclusions 

Our estimates of energy densities in heavy-ion collisions in the CERN SPS 
energy region have been based on the Glauber-type model of tube-on-tube 
interactions. The tube-on-tube prescription of nuclear collisions was comple- 
mented by the introduction of the formation time. 

In spite of our rather conservative choice of model parameters, the energy 
densities we obtained are higher than the critical value known from lattice 
gauge theory. Thus our results support the claim that the threshold for 
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Figure 4: The band in 2La, 2L b plane corresponding to e c = 0.7 ± 
0.35 GeV/fm 3 , calculated for t = 1 fm/c and Ay = 0.5. 

quark-gluon plasma formation has been reached at the CERN SPS |4*H] . 
This statement requires a few comments. 

(i) The lattice results |2*2*ll2*3*] on the critical energy density still have rather 
large error bars. According to [22] e c = (6 ± 2) x (173 ± 8) MeV 4 which 
leads to e c = 0.7 ± 0.35GeV/fm 3 . In addition to statistical errors on 
T c , there are also systematic errors of a comparable size due to extrap- 
olation to the chiral limit. The uncertainty due to statistical errors 
is illustrated in Fig. 0] where we plot the band between the curves for 
e c = 0.35GeV/fm 3 and e = 1.05GeV/fm 3 corresponding to our results 
for Ay = 0.5 and to — 1 fm/c, together with the curves obtained in the 
BO and KLNS models. The curve corresponding to e c = 1.05 GeV/fm 3 
(that means 1.5 times 0.7 GeV/fm 3 ) starts to approach the BO and 
KLNS models. In view of this uncertainty it becomes unclear whether 
or not critical energy density is reached in S+S collisions, but central 
collisions of In+In and heavier systems appear to be on the safe side. 

(ii) We have tacitly assumed that the local energy density as calculated 
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above for a small size system can be compared with lattice results cor- 
responding to an infinite system. 

(iii) The model we are using is based on a phenomenological picture of pro- 
duction of (mostly) soft hadrons. Although the model has been tested 
by comparison with the data on hadron production, the true dynamics 
of the process might be somewhat different. It is for instance possible 
|49j . that the first stage of the collision is dominated by production of 
gluons with momenta of about 0.6-1.0 GeV/c in nucleon-nucleon in- 
teractions and the system — depending on the energy density — either 
hadronizes or approaches kinetic equilibrium. In the latter case the for- 
mation time as used in Eqs. (11) and (12) should rather correspond 
to the thermalization time. This can be roughly estimated as the time 
(in the c.m.s. of colliding nucleons) required for the emission of a few 
softer gluons by the harder ones originally produced. Since the emission 
of a gluon with momentum of kT c takes about 1 fm/c, the approach to 
equilibrium may take 2-3 fm/c. Larger formation time would shift our 
results in Fig. |2] closer to the curves obtained in BO and KLNS models 
from phenomenological analyses of data on anomalous J/ip suppression. 
This is seen in Fig. |3] where we present the energy density obtained in 
our model for values of formation times of 2 and 3 fm/c. Also note that 
lattice results correspond to the system in equilibrium and very little is 
known about J / ip suppression by partonic systems out of equilibrium. 

(iv) The model used is based on the assumption of purely longitudinal dy- 
namics of the nuclear collision — at least in that part of the collision 
which leads to formed hadrons or the system of gluons close to kinetic 
equilibrium. The data on HBT radii and on the transverse momentum 
spectra in heavy-ion collisions at the CERN SPS indicate the presence 
of a rapid onset of the transverse flow [17J E2] . Transverse expansion 
lowers the energy density (with respect to the one calculated in a purely 
longitudinal dynamics) and this would move our curves closer to those 
obtained in the BO and KLNS models. Unfortunately, without a more 
detailed information about the time evolution of the transverse flow it 
is difficult to estimate the effect. 

(v) It has been shown by the Bielefeld group |2Uj and later confirmed by 
Wong and co-workers jjl] that (2S) and (IP) quarkonia may already 
be dissolved below T c . This affects the observed J/ip suppression via 
dissolution of (IP) quarkonia since about 40% of J/ip's in the final state 
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Figure 5: Curves of the critical energy density calculated for formation (equi- 
libration) times to = 0, 1, 2, 3 fm/c, and Ay = 0.5. 

is due to radiative decays of x's. Hence, anomalous J/ip suppression 
starts at energy density lower than e c The discrepancy between the 
lines corresponding to this lower energy density in our model and the 
BO and/or KLNS models shown in Fig. |3] is yet higher than what is 
presented in the Figure. 

To summarize: we have computed energy densities reached in heavy- 
ion collisions at the CERN SPS in a simple model based on the assump- 
tion of longitudinal dynamics in the stage preceding hadronization or kinetic 
equilibration. For formation times of about 1 fm/c the energy density of 
e c = 0.7GeV/fm 3 is reached in collisions of tubes which are shorter than as- 
sumed in phenomenological models of Blaizot and Ollitrault and of Kharzeev, 
Lourengo, Nardi and Satz. The discrepancy is most likely due to combina- 
tion of the three following effects: larger formation or equilibration times 
than usually assumed, true critical energy density larger than 0.7 GeV/fm 3 
and a possible rapid onset of transverse expansion. 
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